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ABSTRACT 



Context. Determining Lyo- properties of distant galaxies is of great interest for various astrophysical studies. 

Aims. We examine how the strength of Lyo- emission can be constrained from broad-band SED fits instead of relying on spectroscopy. 
Methods. We use our SED fitting tool including the effects of nebular emission, considering in particular Lyo- emission as a free 
parameter, and we demonstrate our method with simulations of mock galaxies. Using this tool we analyse a large sample of U, B, V, 
and i dropout galaxies with multi-band photometry. 

Results. We find significant trends of the fraction of galaxies with Lya- emission increasing both with redshift z and towards fainter 
magnitude (at fixed z), and similar trends for the Lyo- equivalent width. Our inferred Lya properties are in good agreement with the 
available spectroscopic observations and other data. 

Conclusions. These results demonstrate that the strength of Lyor emission in distant star-forming galaxies can be inferred quantita- 
tively from broad-band SED fits, at least statistically for sufficiently large samples with a good photometric coverage. 

Key words. Galaxies: starburst - Galaxies: ISM - Galaxies: high-redshift - Ultraviolet: galaxies - Radiative transfer 
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For various reasons observations of the Lya line in high redshift 
galaxies are of great interest. For example, Lya together with 
other observations can provide useful information on the phys- 
ical properties of distant galaxies, such as their outflows, dust 
content, Hi column density, age and others (e . g. IShaole 



2003 



2009; 



Verhamme et alj|2008h lAtek et alJ feOQS; Pe ntericci et al 
Haves et alJ l20ia I201 ll) It mav also be used as a sig- 



X 



nature to distinguish "normal" stellar populations from ex- 
trem ely metal-poor ones, or even to find Population III stars 
(e.g. ISchaere3l2003l: iNagao et alJl2008l) . Understanding the be- 
haviour of Lya also allows us to clarify the overlap between 
different galaxy types such as Ly man break galaxies (L BGs) 
and Lyman alpha emitters (LAEs) (Verhamme et al. 2008) and 
the influence of this line on selection functions dStanwav et alJ 
! I2008t iReddv & Steidelll2009t) . Finally, Lya observations of dis- 
tant galaxies are also being used to study questions of more cos- 
mological nature, such as galaxy clustering and dark energy, and 
to probe cosmic reionisation beyond z > 6 dHamana et alj|2004c 
Mal hotra & Rhoad j2004tlHifl et a l. 2008). 

i 

Generally, the Lyo- properties of distant galaxies have been 
determined from spectroscopic observations (follow-up or blind) 
or via n arrow -band observations targetted to specific redshifts. 
ICookd (I2009I) has shown that LBGs at z ~ 3 dominated by 
Lya absorption or emission can be distinguished based on broad- 
band photometry. Using SE D modeling techniques taking emis - 
sion lines into account (cf. Schaerer& de Barros 2009, 2010), 
we take this one step further demonstrating quantitatively that 
the properties of Lya emission can also be inferred from broad- 
band observations, at least statistically for large samples with 
sufficient photometric bands. This allows us, for example, to de- 
termine trends of Lya with redshift and other parameters, with- 
out resort to spectroscopy. 



Our paper is structured as follows. The observational data 
and the method used for SED modeling are described in Sect. [2] 
The results are presented and discussed in Sect. [3] Sect. [4] sum- 
marises our main conclusions. We adopt a A-CDM cosmological 
model with H =10 km s _1 Mpc" 1 , Q m =0.3 and £2 A =0.7. 

2. Observational data and SED modeling 

Ide Barros et alJ d201 ll) have analysed in depth a large sample of 
z ~ 3-6 dropout selected galaxies using an up-to-date photo- 
metric redshift and SED fitting tool, which treats the effects of 
nebular emission on the SEDs of galaxies. In their homogeneous 
analysis they determine the main physical properties, such as the 
star formation rate (SFR), stellar mass, age, and reddening. They 
assess carefully their uncertainties, and discuss the evolution of 
these properties with redshift. We here extend these simulations 
to perform SED fits with variable strengths of Lya for the same 
set of galaxies, with the aim of examining whether the available 
photometry allows us to distinguish any trends of Lya with red- 
shift and/or other properties. 

2.1. Photometric data and sample selection 

We ha ve used the GOODS-MUSIC catalogue of lSantini et"ail 
d2009h providing photometry in the U, B435, V606, hi6, zssolp, J, 
H, K, bands mostly from the VLT and HST, and the 3.6, 4.5, 5.8, 
and 8.0 fim bands fro m the IRAC camer a onboard Spitzer. Using 
standard criteria as in IStark et al.i (2009) we have then selected 
U, B, V, and i-drop galaxies. To reduce the contamination rate 
(typically ~ 10-20 %) we have only retained the objects whose 
median photometric redshifts agree with the targetted redshift 
range. This leaves us with a sample of 389, 705, 199, and 60 
galaxies with m edian photometric reds hifts of z p h t = 3.3, 3.9, 
4.9, and 6.0. See Ide Barros et all d201 ll) for more details. 
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Fig. 1. Example of an SED fit for B-drop galaxy best fit with strong Lya. Black solid line: best-fit SED including nebular emission. 
Red crosses: flux in the photometric bands derived from the SED including nebular emission. Blue symbols show the observed 
fluxes or 1 cr upper limits. Green dashed lines show the positions/shape of all the availabe photometric bands. Left: Full SED 
showing all the available photometric data. Right: Zoom on the SED showing the spectral range up to 2.4 fim only (restframe UV 
to B aimer break). 
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Fig. 2. Same as Fig.Q]for an object best-fit with a low EW(Lya). 



2.2. SED models 



Our SEP fitting tool, descr i bed in ISchaerer & de Barrosl {2009) 
and ISchaerer & de Barrosl d2010l) . is based on a version of 
the H yperz photometric redshift code of iBolzonella et al.l 
(2000), modi fied t o taki ng into account nebular emission. 
In de Barros et al. I J20T1 we consider a large set of spec- 
tral templates based on the GALAXEV synthesis models of 



iBruzual & Chariot) ([2003), covering different metallicities and 
a wide range of star formation histories. A Salpeter IMF is 
adopted. Nebular emission from continuum processes and nu- 
merous emission lines is added to the s pectra predicted from 
the G ALAXEV models as described in ISchaerer & de Barrosl 
(2009), proportionally to the Lyman continuum photon produc- 
tion. The i ntergalactic m edium (IGM) is treated with the pre- 
scription of lMadaul d 19951) . 

The free parameters of our SED fits are: redshift z, metal- 
licity Z (of stars and gas), star formation history described by 
the timescale t (i.e. SFR oc exp(-f/r)), the age f* defined since 
the onset of star-formation, and a ttenuation Ay described by the 
Calzetti law dCalzetti et alj2000l) . In addition, we here introduce 
a variable Lya strength described by the relative Lya escape 
fraction fj*Le [0, 1], defined by 



L(Lya) = f* a xL B (Lya), 



(1) 



where L B (Lya) is the intrinsic Lya luminosity of the spectral 
template given by its Lyman continuum flux and the case B 
assumption, and L(Lya) is the adopted Lya luminosity for the 
spectral template (before any additional attenuation with the 
Calzetti law, assumed to affect stars and gas in the same man- 
ner). Values f£~ a < 1 therefore describe an additional reduction 
of Lya beyond the attenuation suffered by the UV continuum. 

In practice we compute SED fits for all combinations of z e 
[0,10] in steps of 0.1, Z = (O.O2=Z , 0.004, 0.001), t = (10, 30, 
50, 70, 100, 300, 500, 700, 1000, 300 0, 00) Mvr, 51 age step s 



ge steps 
I. 2000), 



from to the age of the Universe (see Bolzonella et al 
A v e [0,4] mag in steps of 0.1, and f^ a = (0,0.25,0.5,0.75,1.0). 
Minimisation over the entire parameter space yields the best- 
fit parameters as well as other properties such as the predicted 
Lya equivalent width and UV magnitude. To determine the un- 
certainties of the physical parameters, we use Monte-Carlo sim- 
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ulations by running typically 1000 realisations of each object. 
From this we derive the probability distribution function for 
each parameter/quantity, either for each individual object or for 
(sub)samples. 

Examples of best-fit SEDs of two B-drop galaxies where a 
high (low) Lya equivalent width is inferred from the data are 
shown in Figs. [T]© for illustration. The position and shape of 
the availabe photometric bands are also shown. Several issues 
can be seen from these figures. First, we see that 2-3 bands, un- 
affected by Lye, are available to constrain the restframe UV do- 
main, i.e. the UV slope. This is also true for our higher redshift 
samples, where e.g. the JHK bands probe the UV restframe for 
I-drops. Of course, the photometry is not necessary available and 
deep enough for all objects. Second, in general the contribution 
of nebular emission to the broad-band filters can vary strongly 
from case to case. See Fig. Q] showing an object best fit with a 
strong lines and a strong nebular continuum, whereas for the ob- 
ject shown in Fig. [2] nebular emission is negligible. Note, that 
the first object shows a clear excess at 3.6 //m with respect to 
4.5 jum, which is naturally explained by presence of strong lines 
(Ha and others) in the first filter and few lines in the second for 
galaxies in this redshift range (z ~ 3.8-5, cf.lShim etal] l201lt 
Ide Barros et alj|201lh . In this case strong intrinsic Lya emission 
is also expected, and our fitting method yields a median (mean) 
Lya escape fraction /L yff =0.75 (0.57) and a corresponding Lya 
equivalent width EW(Lya)—93 (71) A. Finally, strong Lya emis- 
sion does not necessarily imply a strong excess in the filter en- 
compassing this line with respect to the next redward filter, since 
intrinsically EW(Lya) is not very large, and since the continuum 
flux starts to be reduced by the IGM at z <: 3. 

Typically Lya contributes to ~ 25-30% of the broad-band 
flux for U- to I-drops if we assume EW rest (Lya) = 60 A and 
compare this value to the filter width. However, this contribution 
is higher if we take into account the fact that the continuum flux 
decreases blueward of Lya. Therefore, depending on the source 
redshift and on the exact IGM attenuation, the relative contribu- 
tion of Lya to the photometric signal can significantly exceed 
the above estimate. In any case, when sufficient filters are avail- 
able, the Lya signal can be determined (albeit with considerable 
uncertainty), as we will show below. 

3. Results 

3.1. Evolution of the Lya fraction with redshift and UV 
magnitude 

We have examined the probability distribution function (pdf) of 
the relative Lya escape fraction f^ a of our samples. Overall it 

turns out that is not well constrained for individual objects. 
For each sample, however, the pdf shows two relative maxima 
close to f£f a « and 1, whose relative importance varies be- 
tween the samples. To quantify this behaviour further, regardless 
of the detailed shape of the pdf, we count the number of objects 
with a relative escape fraction above 40%, and define the corre- 
sponding fraction of Lya objects as rL ya = N(f^ a > 0.4)/iV to t- 
This number should provide a simple estimate of the fraction of 
objects with Lya emission among the total dropout samples. 

In Fig. [3] we plot as a function of the absolute UV mag- 
nitude for the samples of z ~ 3, 4, 5, and 6 galaxies. Two main 
results emerge from this figure. First, we find that the fraction 
of objects showing Lya emission increases with redshift (both 
on average and within each magnitude bin). Second, at each red- 
shift we find that Lya emission is more common in galaxies with 
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Fig. 3. Fraction r^ of galaxies with Lya emission (defined by 
fhya > 0-4) as a function of the absolute UV magnitude (at 1500 

A) derived from our U (green line), B (blue), V (black) and I 
(black dotted) drop samples. 



fainter UV magnitudes. Except at the brightest and faintest mag- 
nitudes, where the number of galaxies is relatively small, our 
trends are significant. The same trend is also found when using 
cuts in equivalent width instead of r^ ya , as shown below. 
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Fig. 4. Same as Fig. [3] for galaxies with a large Lya equivalent 
width (EW(Lya)> 50 A) compared to the observed fraction de- 
rived from follow-up spectroscopy of z ~ 4 (blue po ints) and 5 
(black points) LBGs (data for B and V drops by from lStark et al.l 

(boiom . 
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Interestingly our derived Lya fraction shows the same be- 
haviour with redshift and UV magnitude as obtained from and 
hinted at by several other studies, mostly based on spectroscopic 
observations of z ^ 3 galaxies. For example , using photomet- 
ric criteria derived from spectroscopy, Cooke (2009) showed an 
increasing fraction of objects with strong Lya emission among 
faint z ~ 3 LBG s. The mos t direc t comparison can be made 
with the work of IStark et al.l d2010h . who determined the Lya 
fraction of LBGs from spectroscopy of ~ 400 B, V, I, and z 
dropout galaxies. These authors found indeed the same trends 
just described both with UV magnitu de and with reds hift. To 
compare our Lya with that derived by Star k et al.l d2010l) we ap- 
ply the same criterium on the Lya restframe equivalent width 
(EW(Lya)> 50 A) on our samples. The result, shown in Fig. [4] 
is very encouraging, although some discrepancy seems to remain 
at intermediate magnit udes. The s ame holds when we compare 
our data with those of Stark et al. (201 lj). For our U-drop sam- 
ple we find « 10 + 5 % of objects with EW(Lya)> 20 A at 
M \jy < -20.5, co mparable to ~ 20-25 % of LBGs in the sample 
of Shapley et al. (2003) at similar UV magnitudes. Among z ~ 5 
LBGs lDouglas et all (1201 Ol) find ~ 55% with EW(Lya)> 20 A, 
or 22 + 4 % after correction for spectroscopically unconfirmed 
galaxies. Our values for the same EW(Lya) cut are between ~ 

15-60 % for magnitudes brigther than Muv 20. We conclude 

that the Lya trends derived with our method are in reasonable 
agreement with the available spectroscopic data. 

Other data also supports the trend of a higher Lya frac- 
tion at high redshift. For example, a convergence of the lu- 
minosity functions and number densities o f LBGs and LAEs 
has been noticed b y various studies (e.g. Ouchi et al.l 2008: 
Cassata et al.l 1201 ll) . Using such survevs THaves et al.1 (1201 lH : 
Blanc et all d2010l) have shown that, on average, the absolute 
Lya escape fraction^ increases with redshift from z ~ to 6.5 
in UV selected, star-forming galaxies. Our analysis shows that 
this also holds for LBGs. 



3.2. The distribution of EW(Lya) with magnitude and other 
comparisons 

The "recovered" best-fit values of EW(Lya) for all ~ 1300 
galaxies from z ~ 3-5 are plotted in Fig. as a function of 
the UV magnitude. Our models show clearly - at each redshift 
- an absence of UV bright objects with large Lya equivalent 
widths, whereas for fainter objects EW(Lya) occupies a wide 
range. This behaviour is well known, as it has been found in 
virtually al l samples of LAE and LBGs at z ^ 2 (or above 
3) (see e.g. IShaplev et al.1l2003t lAndo et al.ll2006t iTapken et al." 
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2006; Shimasaku et al. 2006; Ouchi et al. 2008; Pentericci et al 
2009t IStark et all 1201 0h T~alfhough iNilsson et all (120091) have 
questioned the statistical robustness of this trend. Our analysis of 
a large sample of LBGs therefore provides an independent con- 
firmation for the deficit of large Lya equivalent widths at bright 
UV continuum luminosities and for an increasing maximum of 
£W(Lya) toward fainter magnitudes found from spectroscopic 
and narrow-band surveys of LBGs and LAEs. The statistical sig- 
nificance of our result has already been shown in Fig. [3] 

The Lya equivalent widths predicted from our SED fits span 
a reasonable range of values, comparable to those observed in 
various LBG samples (cf. e.g. IStarketaO2010t iPentericci et al.l 
120091: fShapley et~aT1l2003h . with most values between and ~ 
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1 The absolute Lya escape fraction is given by fft a X 10 0AAl w , where 



Aiy a is the attenuation at the wavelength of Lya. 



Fig. 5. Predicted best-fit values of the restframe Lya equiva- 
lent width £W(Lya) as a function of Muv for the U, B, and 
V dropout galaxies studied here. Same colour symbols as in the 
previous Figures. To illustrate the typical uncertainties for the 
individual objects, we show errorbars randomly for every fifth 
B-drop. 



60 A (the typical value for constant star-formation with f^~ a « 
1). A good agreeme nt is found between the EW distribution of 
Shapley et al. (2003) and our results for the U-drop sample. 

Finally we have examined all individual objects for which 
spectro scopic measureme nts of Lya where available. From the 
data of IStark et alJ (l2010h we find 36 objects in common with 
our sample. The predicted and observed Lya equivalent widths 
of these objects is shown in Fig. [6] Overall we find an agree- 
ment within a 68% confidence interval for ~ 60% of the objects. 
The main outliers are 3 i-drop galaxies for which our model 
predicts too large equivalent widths (possibly related to uncer- 
tainties with the IGM treatment, cf. below), and several B-drops 
where the model underestimates £W(Lya). Although encourag- 
ing, further comparisons with larger samples may be helpful. 

We conclude that our SED fitting method relying purely on 
broad-band photometry allows us to derive statistically meaning- 
ful trends of Lya with redshift and UV magnitude. Furthermore 
the trends and the predicted Lya strength appears both qualita- 
tively and quantitatively in good/reasonable agreement with the 
existing observations. 



3.3. Test of our method using mock simulations 

The above comparison of our results with 36 galaxies with avail- 
able Lya measurements can be improved by mock simulations 
allowing a more complete assessment of our method. 

To do this we have constructed a set of synthetic galaxy spec- 
tra with known physical properties including EW(Lya), from 
which we then derive the synthetic photometry including typical 
observational errors. The photometry of these mock galaxies is 
then fitted using our SED fitting tool, and the derived EW(Lya) 
compared to the intrinsic (input) one. In practice we have cho- 
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Fig. 6. Comparison of the predicted Lya (rest-frame) equiva- 
lent width with the observed one for the 36 objects i n com - 
mon with spectroscopic measurements of IStark et al.l ((2010). 
Observational errorbars of ~ 25%, typical for the brighter ob- 
jects, are adopted. A logarithmic x-axis is chosen to optimise the 
readability of the plot. Red symbols indicate the i-drops, black 
V-drops, and blue B-drops. For approximately 60% of the ob- 
jects we find an agreement within a 68% confidence interval. 



Fig. 7. Derived median Lya equivalent width (restframe) ver- 
sus input Lya equivalent width for our mock sample of 500 B- 
drop galaxies. The errorbar indicates the median 68% confidence 
limit on the derived EW(Lya); no offset between the derived me- 
dian and the input value is obtained, which demonstrates that our 
method works at least statistically. A small fraction of objects is 
found at higher EW(Lya), outside the plot. A good agreement is 
obtained for all of them. 



sen to simulate a sample of 500 B-drop galaxies. To construct 
their synthetic spectra we randomly draw their physical parame- 
ters, i.e. varying age from to ~ 0.5 Gyr, star-formation histories 
with variable r (taking the same values as for the SED fits, cf. 
above), variable f^ a e [0, 1], A v e [0, 1], and z e [3.8,4.5]. We 
chose the intermediate m etallicity (Z _ 0.004), since our fits are 
very degenerate in Z (cf. lde Barros et al]|201 ll) . Finally we add 
noise to the synthetic photometry taking the median photometric 
uncertainty of the observed B-drop sample in each band. Our ar- 
tificial sample of galaxies should therefore be representative of 
typical B-drop galaxies. 

The mock sample generated in this manner is then fitted us- 
ing the same procedure adopted for the observed LBGs, allow- 
ing us to verify the accuracy of our fitting procedure. The result 
for the Lya equivalent width is shown in Fig. [7] A good agree- 
ment is found between the derived (output) and input equiva- 
lent widths albeit with a relatively large uncertainty. Examing 
the full probability distribution function of the derived equiva- 
lent widths we find (EW(Lya)-EW in ) = 0.+^ A (within 68% 
confidence), which demonstrates that our method is statistically 
able to recover the input value with no systematic offset. As ex- 
pected the dispersion between the input and output equivalent 
widths increases with decreasing input EW(Lya). Note that the 
above dispersion has been obtained for a sample with a relatively 
low median EW m = 20 A. 

One may wonder whether the accuracy of recovering the 
strength of the Lya line depends on UV slope of the galaxy. If 
this was the case the trends found above with UV magnitude 
could be biased, since the o bserved UV slope ap pears to corre- 
late with this quantity (e.g. iBouwens et a"fll2009h . In Fig. [8] we 



show the relative difference ("error") between the output (de- 
rived) and input Lya equivalent width for our mock sample as a 
function of the UV slo pe fi determined in a s tandard fashion from 
the (i776-Z850Lp) color (IBouwens et al.l l2009). For clarity we sep- 
arate galaxies according to their input equivalent width EWj n in 
this plot. We find no trend of the relative "error" with /3 for ob- 
jects with EWj n > 50 A. For objects with faint Lya emission the 
errors become larger, as expected, but show no clear trend with 
UV slope, at least over the range of typical slopes observed at 
z ~ 4 (J3 ~ -3 to -1). We therefore conclude that our multi-band 
SED fitting method is able to recover the strength of Lya emis- 
sion with an uncertainty which is independent of the intrinsic 
UV slope of the galaxy. 

3.4. Possible biases 

Could the main results obtained with our method, i.e. the in- 
creasing fraction of Lye emission towards fainter UV magni- 
tudes and the increase of Lyo- with redshift among LBGs, be due 
to biases and not intrinsic to the galaxy population? We think 
that this is not the case for the following reasons. 

First, our method is not biased by a possible trend of the 
UV slope towards fainter magnitudes, at least for the objects 
with sufficiently strong Lya emission (e.g. EW(Lya)> 50 A). 
Furthermore our photometric method yields the same behaviour 
as spectroscopic measurements, which should be subject to dif- 
ferent biases, and which generally do not reach the same depths. 
To avoid biases due to the difficulty of d etecting faint Lya emis- 
sion from objects faint in the continuum. lStark et al.l d2010) have 
measured the fraction of LBGs at z ~ 4-6 with EW(Lya) above 
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B-drop mock catalogue 




(EW -EW in )/EW in 

Fig. 8. Relative difference between the output (derived) and input 
restframe Lya equivalent width for our mock sample as a func- 
tion of the UV slope /3 determined from the (i776-Z850Lp) color. 
Blue: objects with EW in > 50 A, green: 50 > EW in > 20, black: 
EW in < 20 A. 



an appropriate threshold, and we have applied the same cuts here 
(Fig. |4|. Recent, deeper observations of z ~ 2-3 LBGs by N. 
Reddy (Private communication) also show a persistent increase 
of the Lya fraction towards faint objects once corrected for this 
bias. 

Finally, these trends could also be biased by the LBG selec- 
tion, which can be affected by the strength of Lya emission (cf. 
Stan wavetal.1 120071 120081; iReddv et aljl2008l: iReddv & StekieU 
20091) As pointed out by IStark et alJ (120101) this bias is min- 
imised by simultaneously targetting "adjacent" drop-out popu- 
lations, e.g. U,B,V, and I-drops as done here. In this way objects 
with strong Lya emission scattered out from one selection win- 
dow will be recovered by the selection of the adjacent drop-out. 
Therefore biases related to Lya affecting the colors used for the 
LBG selection should only influence the approximate redshift 
distribution of galaxies, but should not bias the overall sample 
with respect to the strength of Lya. 

We conclude that the derived trends of Lya with magnitude 
and redshift must be an intrinsic property of UV selected star- 
forming galaxies at high redshift. 

3.5. Discussion 

What explains the observed trends of Lya with UV mag- 
nitude and redshift? Based on radiat ion transfer mo d els an d 
on comparisons with obser vations IVerhamme et al.l (2008); 
ISchaerer & Verhammd d2008l) have shown that absorption by 
dust is the main process responsible for the observed diversity 
of Lya line profiles and strengths among LBGs and LAEs. From 
this they suggested that variations of the dust content with galaxy 
mass, UV magnitude, and redshift should explain the increasing 
fraction of Lya emission both with z and with decreasing bright- 
ness of the objects. Findings of a correlation between extinction 



and t he Lya escape fraction (cf. Atek et al1 l2009t lHaves et all 
2010) and the observed decrease of the average UV attenua- 
tion with incr easing redsh i ft furt her support this picture, even 
quantitatively Have s et all d201 ll) . To examine this further, we 
have checked whether our SED models reveal any correlation 
between ff~ a and Ay (or other physical properties). Beyond the 
one shown above with M\jy, we have not found clear additional 
correlations, the uncertainties/degeneracies being too large. We 
can speculate that future studies of large galaxy samples with 
additional observational constraints, e.g. from deep surveys with 
large sets of filters, should be able to detect such correlations. 

Can our method be used to constrain cosmic reionisation? In 
principle yes, since as shown above, we are able to determine 
the strength of Lya in a similar way as spectroscopic observa- 
tions, at least for a large sample of galaxies. From this we are 
then in principle able to determine the fraction of Lya emitters, 
the Lya luminosity function, or other related quantities, which 
can in turn be used to probe the ev olution of the neutral hy- 
drogen fraction with redshift (see e.g. Malhotra & Rhoads 2004; 
Diik stra et al.l l2007T). However, we should be aware of the fol- 
lowing possible limitations. At high-z (z » 5) the IGM signif- 
icantly alters the SED shortward of Lya. In order to correctly 
predict the Lya strength from broad-band SED fits one there- 
fore also needs the most accurate treatment of the Lya forest 
(in principle of stochastic nature) possible. At the same time 
there are some degeneracies between Lya strength, IGM, and 
also redshift, which can all modify the flux in the broad-band 
filter encompassing Lya. It is therefore clear that our method 
should become less accurate at very high redshift. Future appli- 
cations to larger samples of high-z galaxies observed in as many 
photometric bands as possible, and careful examination of the 
above mentioned degeneracies will show whether our SED fit- 
ting method can provide useful constraints also on Lya during 
the epoch of reionisation. 

4. Summary and conclusions 

Using an updated version of t he Hyperz photometric redshift 
code of iBolzonella et akl (1200 0) addin g nebular emi s sion to the 
spectral templates (see ISchaerer & de Barrosll2009L 120101) . we 
have analysed a large sample of Lyman-break selected galaxies 
at z ~ 3-6 in the GOODS-S field, for which deep multi-band 
photometry from the U-band to 8 fim is available. To our exten- 
sive exploration of the parameter space covering redshift, metal- 
licity, star-formation hist ory, age, and attenuati on whose results 
are discussed in detail in lde Barros et alJ d201 lb . we have added 
a variable Lya strength, described by a relative Lya escape frac- 
tion fff a . We show that significant trends of Lya strength with 
redshift and with UV magnitude can be inferred from broad- 
band photometric observations of large samples of galaxies us- 
ing our method. The validity of our method is also demonstrated 
by our simulations of mock galaxies. Our method relies on fitting 
all the photometric bands constraining the UV and optical rest- 
frame SED. However, even with just 3 band s some Lya tre nds 
may be detectable, as e.g. shown already by ICookel d2009l) for 
z ~ 3 galaxies. 

We find an increase of the fraction of LBGs showing Lya 
emission both with redshift, and towards fainter UV-restframe 
magnitudes at each redshift. Furthermore, we also infer a de- 
crease of the average and the maximum Lya equivalent width 
toward bright objects. Our results are fully compatible with 
the trends found from sp ectroscopic observations of individ- 
ual LBGs and LAEs (e.g. lAndo et"aT]l2006t lOuchi et alJl2008L 



Daniel Schaerer, Stephane de Barros, Daniel P. Stark: On Lya emission in z ~ 3-6 UV-selected galaxies 



7 



Star k et al.ll2010l) . and with o ther more global tr ends of the Lya 
escape fraction with redshift (lHaves et al.ll201 ll) . 

In principle our method should also be applicable to derive 
the Lya properties of z > 6 galaxies, of particular interest for 
studies of the history of cosmic reionisation, and in others con- 
texts. At least we have shown that even broad-band photometric 
observations can reveal some information about Lya emission in 
high redshift galaxies. 
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